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The equation of state of the string model with linear strings comes close to describing the lattice 
QCD results and allows for the E /N ~ BTq = 1 GeV relation found in phenomenological statistical 
model. The E/N value is derived from the zero pressure condition in quark matter and is a fairly 
general result. The baryochemical potential dependence of the hadron gas can be met if it is re- 
- - - , interpreted in the framework of an additive quark model. The conclusion is reached that stringy 

9^ ■ models explain the E/N — 6To relation naturally and independently of the value of the string 

' tension. 
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The phenomenological success of the statistical hadronization model has been emphasized, questioned and criticised 
Oh repeatedly over the years. All experimental results of heavy-ion collisions on particle yields at energies ranging from SIS 
^ to RHIC are consistent with results falling in a narrow stripe in the parameter space of temperature and baryochemical 
rtj ' potential. There are several descriptions of this so called chemical freeze-out curve, for recent papers on this topic see 

^: mm. 

In this paper we intend to explain the rather high value of the energy per particle, E/N « 1 GeV for a system 
' of quarks and giuons with an equation of state containing a term depending on the color density. As a matter of 
^ ' fact the convex shape of this curve (cf. Fig|3]) is in agreement with several thermodynamic approaches incorporating 
00 ' fermions and bosons as ideal gases. The main puzzle is the quantitative value, since 1 GeV is six times higher than 
0^ ■ the associated hadronization temperature, Tq « 167 MeV. We aim to understand, how this value can be obtained 
0^ ■ starting from a massless quark-gluon plasma by adding a color-density term. 

, On the hadronic side, considering massive matter, elementary non-relativistic thermodynamics leads to E/N — 
T— I ' m ~\- 3T/2. With the known values cited above this requires an average mass in the range of the rho meson mass, 

' m ~ 750 MeV. This value satisfies m ^ T, so the non-relativistic approximation turns out to be acceptable. 
25 However, such a high mass on the quark matter side cannot be a constituent mass, it can only be the result of 

strong interactions. A well known example for treating this interaction as a mean field, the original MIT bag model, 
with p = kT^ ~ B, e ^ 3kT'^ + B and n « kT^ pressure, energy density and particle density, respectively, would 
allow for E/N = e/n < AT only, upon p > 0. Moreover the equation of state of the model quark matter has to 
be in accordance to results obtained from lattice QGD, the only model independent non-perturbative approach from 
field theory to the equation of state at present. We present a simple thermodynamical model of massless quarks and 
giuons whose interactions are described by a free energy contribution motivated by strings. This stringy matter will 
be studied at vanishing pressure for the hadronization curve and at high temperature for the lattice equation of state. 
It is an interesting conclusion that such stringy models do explain the E/N = 6Tq relation easily and, remarkably, 
independent of the value of the string tension. 
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II. STRINGY THERMODYNAMICS 



There are several ways to treat corrections to an ideal gas equation of state (eos) . Each model identifies a physical 
picture in which the leading order interactions are calculated. The lattice QCD and pure Yang Mills eos approach an 
effective massless ideal gas limit at high temperatures, but it deviates vastly from it around the color deconfinement 
temperature. In this paper we test the string model picture for an interacting gluonic plasma, introduced in 

We assume that a however decreasing proportion of color charges are still connected by strings at high temperature, 
above the crossover value to deconfinement T > Tc. This assumption is supported by some simple and general 
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properties of the thermal distribution of relative momentum squared between pairs of massless particles 
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Pm = (-5 (g^ - 2E,E,{1 - cos^))) = '-^^ . (1) 



This distribution has a non-perturbative contribution to the equation of state stemming from relative momenta below 
a fixed value, (Q^ < with say Awl GeV): 
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F = j P{Q')dQ\ (2) 



This quantity, not having another energy scale than the temperature in case of massless particles is given by a scaled 
integral 

F= j fix)dx. (3) 



For high temperatures T ^ A this proportion is approximately given by F « /(O) A^/T^, for low temperature, T <C A 
by _F w 1 due to normalization. Since 

T \ ' 



is non-zero (cf. eq. [T|), there is always a non-perturbative contribution to the equation of state at any temperature, 
it just becomes subleading order in the pressure. Nevertheless in the interaction measure, (e — 3p), such an O(T^A^) 
term mixes to the leading order. 

A. Equation of state with strings 

A general class of eos is obtained by taking into account a term in the free energy proportional to a fractional power 
of the density. Here we consider a " color density" , a weighted sum of the number of particles: 



^mci. (5) 



Colorless objects do not pull strings, nor do they take part in screening (ending) them. Furthermore different color 
charges like quarks or gluons may have different effective string constants this effect is taken into account in the 
factors Ci. 

The contribution to strings of an average length of 

{£) cx c-'^ (6) 

shall be multiplied by the density, with 7 being a fractional power between zero and one. For straight strings in three 
dimension 7 = 1/3. 

Since we consider density dependent modifications of the eos at a given temperature, we shall use the free energy 
density, denoted /, as the fundamental potential. We consider 

/ = /(n.,T) + -^ci-^. (7) 

1-7 

Here the coefficient A comprises eventual average geometrical shape factors besides the string tension, so it cannot 
be taken from the Regge slope of meson resonances directly. We shall rather fit it to the lattice eos later. The / is a 
general free energy density, for describing a QGP to be specified later. 

The chemical potentials associated to the component i are given by the derivatives with respect to n^: 

tJ-i ^ fj-i + Ac"'^ a (8) 
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where fli = df/dui follows from /. We use so far the additive form of eq.([7]). The chemical equilibrium establishes, 
if possible, at constant values of this chemical potentials. These values are determined by the conserved charges of 
the components. For the sake of simplicity we consider here the baryon charge only (zero for gluons, 1/3 and —1/3 
for quarks and antiquarks, respectively), but further quantities may be introduced into this scheme easily. From the 
equations 

Mi = (9) 

the corresponding number densities can be expressed 

71°^ = Vi[T,iiB]qi,Ci). (10) 

In this chemical equilibrium situation the color density defined in eq. ^ takes its equilibrium value expressed by the 
Vi-a. 

While the real solution of this system of equations can be involved in the general case, in some particular systems, 
e.g. for a massless ideal gas with strings in the Boltzmann approximation, it can be given in analytic form (see later). 
A remarkable general property, however, can be derived without the explicit form of this solution. This property 
related to the hadronization of the QGP is the energy per particle at the edge of the stability: when the pressure 
vanishes. For this purpose we obtain the entropy density, 

s--% = s. (11) 



the pressure 



and the energy density 



P = V/i.n.-/=p--^ci-^ (12) 
^ — ' — 7 



e = f + Ts = i+-^c^-\ (13) 
1-7 



At the endpoint of the mechanical stability p = po = Q and therefore 

7^ .1-7 



P^ = -^c'-\ (14) 
1-7 



From this the color weighted density can be obtained at this point 

L -' 
A7 



Co = ( —J-^P^ ) ' (15) 



and the energy density is expressed as 



eo = eo+Po/7- (16) 



This is a remarkable result. For massless constituents only p — T'^(p{fiB /T) and therefore e — 3p. For straight strings 
7 = 1/3, so we arrive at eo — 6po- As long as the Boltzmann approximation is applicable, for an ideal gas mixture 
p — nT with n — Ui and one concludes that 

I = eo/n = 6T0. (17) 

With To — 167 MeV fitted to hadronization data (and predicted by lattice QCD as the crossover temperature) one 
obtains E /N = 1 GeV at this point. This derivation applies to the quark-gluon side of the hadronization curve. 

B. Chemical equilibrium with strings 

From eqs.(l8|) and ([9]) follows that the equilibrium number densities in general satisfy 

= "-i,A=o(^' lii^B - Ac^^a). (18) 
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Summing with the color weight factors, Ci we arrive at an implicit equation for the effective color (string source) 
density: 

c = ^Ciri!'^j^^„{T,qifiB ~ Ac'"^Ci). (19) 

Its solution in general depends on the temperature T, baryochemical potential /is and the specific color and baryon 
charges Ci and qt. 

Since this implicit equation contains the unknown, c, on its right hand side as a given correction to the chemical 
potential only, we are able to deduce that there is no solution beyond a critical point. This point is characterized by 
the fact that the derivative of eq. p9|) with respect to c is also satisfied, 

1 = jAc~^-' ^(^' - Ac-^c,). (20) 

We denote the color density in this endpoint where the chemical equilibrium ceases to be possible by c^;. A further 
question is how this endpoint of possible chemical equilibrium solutions is related to the zero pressure boundary. 

In the particular case of Boltzmann-like dependence on the chemical potential dui/dfj, = rii/T (i.e. the exponential 
of /i/T is a factor in the number density). We note that for a Boltzmann system with all c^'s equal to either one or 
zero = Ci and therefore the ratio of eqs. ipO)) and gives a simple condition for the critical color density c = ce- 

7^c^^ = l. (21) 
In the same approximation one can obtain an analytic solution to the chemical equilibrium problem. We get 

c=(j2^^(T)e''''"/^^ e-^^-y^ (22) 

with 

X.(T) = <\^omO). (23) 

The sum in the brackets is a c-indepcndent factor, we denote it by (p(T, /is) = T^a^{fiB,T). The solution of eg. ([22]) 
is given by 

.= ^ = -m-^.-) (24) 

with W{x) being Lambert's W-function defining the solution of the transcendental equation We^ — x. The form 
occurring in eq. ((24|) . —W{-'x) has a real solution for positive x values only if x > 1/e. The argument x — is the 
endpoint of stringy thermodynamics, the endpoint temperature is connected to the effective string tension A as 

jA^{TE,fiB)-' = e-\ (25) 

IE 

For considering massless particles as sources and endings of strings the factor (p — a^T'^ has a special form: a{fiB,T) 
can be a function of the ratio (J-b/T only. At vanishing baryochemical potential, fiB = 0, a is a constant. In this case 
ea. p5|) can be resolved for the endpoint temperature: 

Te = {-fAea-^'')^ . (26) 



In particular for 7 = 1/3 it is Te — \J Ae/3a and a — a{fj,B/TE). This implicitely determines the endlinc ij.b{Te)- 



III. MASSLESS BOLTZMANN GAS WITH STRINGS AT ^iu = 

In the foUowings we investigate the one component massless Boltzmann gas at vanishing net baryon density. In this 
case a is constant and the pressure and energy density are expressed easily with help of explicit chemical equilibrium 
solution 
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with 

/rp \ 37+1 

^ = (^^ j . (28) 

Since for ideal Boltzmann gases the Boyle-Mariotte law holds, p = cT, and for vanishing mass e = ip the total 
pressure is reconstructed as (cf. ea lT2|) : 



p^cT[l-^] (29) 



where z is taken from the solution (j24p . The energy density is given by 



^ = ^^(3 + :^(r3^) (30) 



A. High temperature asymptotics 

When comparing with lattice results the high temperature asymptotics of the stringy massless Boltzmann eos 
described by the eqs. (|24l29l30p implicitly is particularly important. In our model T ~> 00 is realized by the x ^ 
limit (we consider 7 G (0, 1) only). The Lambert W- function has the series expansion z = —W{~x) — x + + . . . 
and hence the high temperature (low-x) color density behaves like 

c^a^T^ (l--a; + ...) , (31) 



the pressure like 



and the energy density like 



The interaction measure is given by 
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e = a^r^ ( 3 + 4t — \x + ...). (33) 
7(1-7) 



e - 3p = a^r^ ( 4t-^2; + ...). (34) 
7(1 - 7) 



This analysis reveals that the high temperature asymptotics of the pressure is T'* as it should be; this fixes the value 
of a? by comparing to lattice eos. (5.21 for 2+1 flavor QCD, 1.7 for pure SU(3)). 

The power 7 can be fixed from the subleading terms (and leading term in the interaction measure), xT^. Using 
eg. (1^5)) such terms scale like T^-^i . In order to agree with lattice eos results, i.e. xT'^ ~ the only possible choice 
we are left with is 7 = 1/3. This agrees with our statistical argument for leading non-perturbative effects on the eos 
at T > Tc, and in particular favors the straight string geometry in quark matter. 



B. The zero pressure point as reference 



Since on the lattice only positive pressure states can be achieved by Monte Carlo techniques, the zero pressure 
point. To, is a more physical reference than the endpoint of chemical equilibrium. This occurs at = 1 — 7 = 2/3 and 
due to z = xe"" at xo = fe~^/^ « 0.34. On the other hand xq = e~^{TE/ToY due to eq.(l28]), so we arrive at Tq/Te = 
(xoe)^^/^ « 1.04. In the stringy massless Boltzmann quark matter approximation these two temperatures are quite 
close to each other. Reference to the zero pressure point can be taken by substituting x — xo{To/T)'^ w 0.34(To/r)^ 
in general. 
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FIG. 1; Equation of state from lattice and from the massless Boltzmann stringy model. Pure SU(3) gauge theory results 7] to 
the top, 2+1 flavor QCD results [1,1, [TO] 

to the bottom. The scaled pressure (left sides) and interaction measure (right sides) 

curves are compared. 



With respect to Tq the scaled pressure and interaction measure are given as 



_P_ 

-3p 



1 - 



3z 

y 



(35) 



This reveals a factor of 9xo ~ 3.04 between the high-temperature asymptotic values of the scaled interaction measure 
(e — 3p)/T^Tq and scaled pressure p/T'^- Although this gives a qualitatively correct picture of the quark matter eos, 
unfortunately for a quantitative comparison Monte Carlo data are still spread over a remarkable range (cf. Fig[T]). 

The whole temperature dependence of the stringy eos in Boltzmann approximation works surprisingly well for the 
pressure, but is only qualitatively correct for the interaction measure. This we plot in Fig[TJ The lattice data show 
the most sudden jump in the T^-scaled interaction measure to a nearly constant value. For comparison the stringy 
model has a somewhat milder rise. Quantitatively the constant value of the scaled interaction measure differs from 
the lattice result in the pure SU(3) case [2] as well as probably also in the full QCD case [1, [§, The massless 
Boltzmann approximation is clearly an oversimplification. 
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C. Degenerate Fermi gas at T=0 



Another analytically solvable case is represented by the degenerate Fermi gas at zero temperature with stringy 
contributions to thermodynamics. In this case, considering only one type of fermions with color charge = 1 the 
chemical potential fixes the Fermi energy. Let us consider the case of massless fermions, than the Fermi momentum 
is also n and the chemical equilibrium relates the color density to this potential via 

c=^{,-Ac-y^)\ (36) 

with degeneracy factor d — (25 + l)NcNf — 12 for light quark matter. 

This equation is analytically solvable. Introducing the parameter B = (Btt^/c?)^/'^ we have a second order algebraic 
equation for the variable c^^^: 

Bci/3 + Ac-i/3^/i. (37) 

The stable solution belongs to the higher value of c (because if the Ihs of eq. ((37|l exceeds fi then chemical processes 
tend to decrease c): 



It is easy to see that the critical value of the chemical potential is given by 

fiE - 2VAB, (39) 

for fj, < fiE there is no equilibrium, the strings pull the matter infinitely. 

For obtaining the pressure and energy density we use the corresponding expressions for a degenerate Fermi gas: 

4 

^ I * — 1 \ 

P = Mr 

3p (40) 



e — 3i 

Upon using eg. ([55)1 this leads to the total pressure 



i^c- -v4c2/3, (41) 



and to the energy density 

3 3 
e = —I1C+ — 



Ac^^\ (42) 



At zero pressure, p = 0, the energy per particle becomes e/n = fj, and the value of the chemical potential at this point 
can be obtained from eq. ([4T|) as satisfying ^qcq = 3Acq^^. Together with eq. ([36|) this leads to 



Mo = 3^AB/2 = « IMfiE (43) 

at the color density Cq = {2A/B)^/^. If this meets the E/N = 1 GeV fine, then at T = = 1 GeV. 

Since the endpoint and the zero pressure points are close to each other both at zero baryochemical potential and 
at zero temperature, we conjecture, that this is the case all over the hadronization line. The baryochemical potential 
at the endpoint of the stringy QGP is given by the minimum value of the curve 

SAB ^, 

fiB = -, — + Skp (44) 

with kp being the Fermi momentum (for massless fermions it is also the Fermi energy, i.e. the chemical potential 
of the quarks). Its minimum is located at fcp = VAB, the value being /i^ = 6\/ AB. However, regarding the 
hadrons fitted in the Statistical Model as objects made in the constituent (additive) quark model, one has to consider 
fig = 3/co = 3^/ AB. This value is exactly the half of the previous one. Fig[2] plots this relation. 
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FIG. 2: The baryochemical potential as a function of the Fermi momentum of quarks at T = in the stringy QGP. The 
minimum point represents the endpoint of the stringy model, the corresponding additive quark model value is the half of it. 
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FIG. 3: The zero pressure line of the stringy, massless QGP (labeled hy q — l/3)with 2 and 3 quark flavors. For comparison 
the Statistical Model results (boxes), fit by T = 0.167 — 0.l39fi% ~ 0.053/x^, Ref.jj] are indicated. We also plotted the stringy 
QGP hadronization lines with halved baryochemical potentials roughly corresponding to an additive quark model interpretation 
of hadrons (labeled by g = 2/3). 



The corresponding hadronization curves are plotted in Fig[3) The p — line (which is close to the endline of 
the stringy thermodynamics within a few per cent) is indicated for a quark gluon plasma with massless gluons and 
Nf = 2 and Nf = 3 massless quark and antiquark flavors, respectively. For comparison the Statistical Model results 
are indicated by filled boxes. Following our previous comment about the additive quark model interpretation of the 
stringy QGP results, the same hadronization lines are also plotted as functions of three times the quark Fermi level 
instead of the baryochemical potential. In this second case the coverage with the hadronic fit results is intriguing. In 
Fig|4]wc plot the energy per particle along the hadronization line; it is roughly constant for moderately low values of 
the baryochemical potential, but then decreases somewhat. 

Summarizing the stringy eos with 7=1/3 comes close to describing the T > Tq lattice QCD results and at the 
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FIG. 4: The energy per particle for a stringy, massless QGP with 3 light flavors. Until the value fiB « 0.8 GeV the E/N ratio 
does not change more than 10 per cent. 



same time offers an explanation for the E/N « 6To = 1 GeV hadronization condition found in phenomenological 
statistical model. While the massless Boltzmann approximation should be dropped in the view of lattice eos data, 
the E/N value at zero pressure is a much more general result; it hopefully survives as a possible physical picture 
for the quark matter side at hadronization in rclativistic heavy ion collisions. Admittedly there is no reason for the 
E/N ratios being equal in the hadronic and quark matter, but it is not easy to imagine, how this value could be 
doubled from 3To to 6To at the hadronization. The baryochemical potential dependence of the hadronic side fit by 
the Statistical Model can be met if the chemical potential is re-interpreted in the framework of an additive quark 
model (parameterized by the half of the stringy QGP value). 
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